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Abstract 
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^ | Elastic nd scattering is considered within the Alt-Grassberger-Sandhas 

(AGS) formalism for various nN input data. A three-body resonant state 
is found close to the nd threshold. This resonance is sustained for different 
choices of the two-body nN scattering length a V N- The position of the res- 
onance moves towards the nd threshold when Rea v ]y is increased, and turns 
J> ' into a quasi-bound state at Kea n N = 0.733 fm. 

in' 
m 

O ■ PACS numbers: 25.80.-e, 21.45. +v, 25.10,+s 

One of the main questions of r^-meson physics concerns the existence of ^-nuclei, a 
Q\ '. possibility predicted in a pioneering work by Haider and Liu [|IJ. In this context it is worth 
mentioning that according to |2j the mean free path of r\ mesons in a nuclear medium is about 
2fm, i.e., less than the size of a typical nucleus. A necessary condition for the existence of 
?7-nuclei, hence, appears to be satisfied. A further indication is the observation of pionic 
nuclei (deeply bound pion states inside nuclei) 0, which were predicted theoretically in j|. 

Recently, in a number of theoretical investigations based on quite different approaches, 
such as the mean field approximation ||, the optical model 0, and few-body calculations 
|7]||, the idea of ^-nuclei was given a rather firm ground. Experimentally, after the discour- 
aging attempt of Ref. ||, certain evidence for the existence of //-nuclei was noticed by two 



groups [|Hjyill . An additional indication is given by the enhancement of near-threshold r\ 
production in the reaction np — > dr] as reported by the Uppsala group [Q . This observation 
is most likely to be understood as the effect of a near-threshold bound or resonant rjd state. 



Moreover, it was suggested [T3] that even quite exotic systems containing an rj meson and a 
hyper on (?7-hypernuclei) may exist. 

Being of interest by itself, the existence of 77-nuclei would also shed new light on various 
fundamental problems of particle physics. For example, the study of 77-nuclei would give 
a clue for understanding the possible restoration of chiral symmetry in a nuclear medium, 
or its partial restoration which may occur at normal densities ~ 0.17fm -3 [|I4"|-|T6[| . A 
further aspect, which makes such systems interesting for other fields of nuclear physics, 
is the modification of the two-body r)N force by the nuclear medium. Being enhanced 
due to the resonant character of the rjN interaction, general properties of such modifications 
should become particularly evident in this case. Especially, the structure of the Sn resonance 
embedded in a nuclear medium may be studied. This would exhibit certain details of the 



effective Lagrangian models of such a resonance [17 , and shed some light, for example, on 
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chiral models which suggest a reduction of its partial width for decaying into the r/N channel 
Wj , or assume Sn to be a quasi-bound state of the K meson and a £ hyperon fl8fl . The 



interesting suggestion of considering Sn as a joint manifestation of threshold (casp) and 
resonance phenomena fl9j could also be checked. 



In a sense, the problem of choosing the most adequate model for describing the Sn 
resonance is similar to the usual difficulty of choosing the "correct" two-body potential. A 
reliable judgment on such a choice can be given when this potential (or the resonance) is used 
in a few-body system where its off-shell properties are exposed. This requires a treatment 
based on rigorous few-body theory. The Alt-Grassberger-Sandhas (AGS) equations PU| . 
which are employed in this work to calculate the //-deuteron elastic scattering amplitude, 
belong to this category. 

The Faddeev-type coupling of these equations guarantees uniqueness of their solutions. 
Moreover, as equations for the elastic and rearrangement operators U aa and Up a , respec- 
tively, they are well-defined in momentum space, providing thus the desired scattering am- 
plitudes in a most direct and technically reliable manner. The advantage of working with 
coupled equations is not only suggested by questions of uniqueness, but also by the rele- 
vance of rescattering effects. Indeed, in recent calculations of //-photoproduction from the 
deuteron it was found that rescattering terms give a significant contribution to the corre- 



sponding amplitude [21 



Let the r\ meson be denoted as particle 1 and the two nucleons as particles 2 and 3. 
With the momentum qi of the meson relative to the center of mass of the nucleons in the 
deuteron state \ipd), the free rjd channel state is given by 

l^;qi) = l^)|qi) , (1) 

its normalization being chosen as 

(V^q'^qi) =%'i-qi) • (2) 

In terms of the AGS transition operator Un the 7]d elastic scattering amplitude is represented 

as 

/(q'i,qi) = -(27r) 2 M 1 (^;q' 1 |f/ii(^)|^;qi) , (3) 

with the on-energy-shell conditions Iq'J = |qi| and z = + qf/2Mi + iO. Here Mi is the 
reduced mass of particle 1 and the (23) subsystem, 1/Mi = 1/mi + l/(m2 + m.3). 



The elastic transition operator Un satisfies the set of AGS equations 

3 

U Pa {z) = (1 - 6pa)Go \z) + £(1 - 6^)T 7 (z)G (z)U^ a (z) , (4) 

7=1 

with Gq(z) being the free resolvent (Green's operator) of the three particles involved. This 
set of equations couples all 3 x 3 elastic and rearrangement operators U aa and Up a . Here 
each of the subscripts runs through the values 1,2 and 3, indicating the two-fragment par- 
titions (1,23), (2,31) and (3,12) respectively. Therefore, Un describes the elastic transition 
1(23) — ^1 (23), while U21 represents the rearrangement process 1(23)— >2(13). 
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In accordance with Eq.([|), we denote by q a the momentum of particle a relative to 
the center of mass of the (^7) subsystem, and by M a the corresponding reduced mass. The 
internal momentum of this subsystem is denoted as p a . Conventionally, this complementary 
notation is also used to label the two-body T-operator t a (z) of the (^7) pair, for instance 
t\{z) = tjvAr(^)- It should be noticed, however, that it is not this genuine two-body operator 
which enters the AGS equations, but the operator 

T a (z)=t a (z-q 2 j2M a ) , (5) 

which is to be understood as the two-body operator embedded in the three-body space, 
with the relative kinetic energy operator q 2 x /2M a of particle a being subtracted from the 
total energy variable z. Considered in momentum space, Eq.@ thus reads 

(p' a ,q' a \T a (z)\p a ,q a ) = 5(q' a -q a )(p' a \t a (z-q 2 j2M a )\p a ) . (6) 

Since we are interested in rjd collision the subscript a is fixed to 1. Instead of the 9 
equations of system (^) we therefore have to consider only the three equations 

U n (z) = T 2 (z)G (z)U 21 (z) +T 3 (z)G Q (z)U 31 (z) , 

U 2l (z) = G \z) + T 1 {z)G {z)U 11 {z) + T 3 (z)G (z)U 31 (z) , (7) 

U 31 (z) = G \z) +T 1 (z)G (z)Uu(z) + T 2 (z)G (z)U 21 (z) , 

involving the operator U\\ which determines the elastic amplitude (0). In momentum space 
the system (0) consists, after partial wave decomposition, of three coupled two-dimensional 
integral equations. It is customary to reduce the dimension of these equations by approxi- 
mating or representing the two-body T-operators in separable form, 

t a (z) = \Xa)T a (z)(x a \ , « = 1, 2, 3 , (8) 

or, according to (|5]), by 

T a {z) = \ X aK{z-ql/2M a ){xa\ , a = 1,2,3. (9) 

Inserting this representation in the AGS equations (Q) and sandwiching them between 
and \xa) they take the form 

3 "2 

X Pa {z) = (1 - 5 0a )( X p\G Q (z)\ Xa ) + XX 1 - W<X/j|Go(*)|X7>t 7 (* - • (10) 

7=1 ZiW 7 

The three-body operators Up a , hence, are replaced by the effective two-body operators 

X 0a (z) = ( X/ 3\G (z)Up a (z)G {z)\ X a) (11) 

which act exclusively on the relative momentum states |q a ). After partial wave decomposi- 
tion we then end up with one-dimensional integral equations. 

The form-factor \xi) in the iViV operator t\ is related to the deuteron wave function 
according to 
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IVtf Qi) = G (g?/2M a + £ d ) |xi)|qi) , (12) 

where is the deuteron energy. The on-energy-shell matrix element in (13), thus, can be 
written in the form 

/(q'i,qi) = -(27r) 2 M 1 (q' 1 |X 11 (^)|qi) . (13) 

In other words, the set of equations fllOD provides the elastic amplitude we are interested in. 

In the present case there are further simplifications caused by the identity of the nucleons. 
Indeed, the momentum representation of X 3i , T3, |xs) and of X 2 i, t 2 , \x 2 ) respectively are 
of the same functional form. This reduces the three equations involving X\\ to the following 
pair 

X u (z) = 2(xi\Go(z)\x 2 )T2 (z - q\/2M 2 ) X 21 (z) , 

X 21 (z) = (X2\G {z)\ X i) + (X2\G {z)\xi)r 1 (z - qj/^h) X 11 {z) (14) 
+ (X2\G (z)\x 2 }T2 (z - c£/2M 2 ) X 21 (z) . 



The S'-wave nucleon-nucleon separable potential is adopted from Ref. |22[ with its pa- 
rameters slightly modified to be consistent with more recent NN data (see Ref. ||). The 
77-nucleon T-matrix is taken in the form 



t vN (p\ P ; z) = (/ + a 2 )- 1 ? *. r/9 A p 2 + ^r 1 (15) 

(z - E + iT/2) 

consisting of two vertex functions and the Sn-propagator in between |7|. It corresponds 
to the process r]N — > Sn — > rjN which at low energies is dominant. The range parameter 
a = 3.316 fm _1 was determined in Ref. PB[ , while Eq and T are the parameters of the Sn 
resonance pl[| , 

E = 1535 MeV - {m N + m v ) , T = 150 MeV . 
The strength parameter A is chosen to reproduce the 7/-nucleon scattering length a V N, 

a\E - iT/2) 

It is customary to use complex a v N the imaginary part of which accounts for the flux losses 
into the irN channel. The value of a V N is not accurately known. Different analyses p5 
provided for a V N values in the range 

0.27 fm < Re a vN < 0.98 fm , 0.19 fm < Im a vN < 0.37 fm . (17) 

Recently, however, most of the authors agreed that Ima^jy is around 0.3 fm. But for Rea^ 
the estimates are still very different (compare, for example, Refs. [^6| and ||27|| ). We, there- 
fore, fixed Im a n N to 0.3 fm and did calculations for several values of Re a V N within the above 
interval. 

We solved Eqs. (|i~4D for rjd collision energies varying from zero (r/rf-threshold, z = E^) up 
to 22 MeV. As is well known (see, for example, f2"8fl), the kernels of Eqs.(|l^), when expressed 
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in momentum representation, have logarithmic singularities for z > 0. These singularities 
stem from the Green's functions which in their denominators involve two terms g 2 /2M Q and 
corresponding to the pair of Jacobi momenta. The singularities appear after angular 
integration, and their position depends on the values of p a and q a (the so-called moving 
singularities). In the numerical procedure, we handle it with the method suggested in Ref. 
|29fl . The main idea of this method consists in expanding the unknown solutions (in the area 
covering the singular points) in certain polynomials and subsequent analytic integration of 
the singular part of the kernels. 

The results of our calculations are presented in Figs. |T]-[^ and in Table |. In Figs. [I] 
and ^ the energy dependence of the r/d phase-shifts for five different choices of Re is 
shown, namely, for 0.55 fm, 0.65 fm, 0.725 fm, 0.75 fm, and 0.85 fm. The larger this value, 
the stronger is the r/N attraction. The change in the character of these curves, hence, reflects 
the growth of the attractive force between the n meson and the nucleon. The lower three 
curves for Re 8 n d corresponding to the smaller values of Re start from zero, the two curves 
corresponding to the strong attraction start from ir. According to Levinson's theorem, the 
phase shift at threshold energy is equal to the number of bound states n times 7r. We 
found that the transition from the lower family of the curves to the upper one happens at 
Re 0^=0.733 fm. Therefore, the 7/N force, which generates Rea^ > 0.733 fm, is sufficiently 
attractive to bind r\ inside the deuteron. This is the first conclusion of our calculations. 

The second conclusion concerns the peaks in the energy dependence of the total elastic 
cross-section (see Fig. f|), indicating that a resonance appears in the rjd-system. Of course, 
not every maximum of the cross-section is a resonance, but the Argand plots, shown in 
Fig. |5], prove that the maxima we found are resonances. Their positions and widths for 
various choices of Rea^ are given in Table |. It should be noted that, while the resonance 
energy is determined in our calculations exactly (as the maximum of the function sin 2 Re 8 v d), 
the corresponding width is obtained by fitting the cross-section with a Breit-Wigner curve. 
Therefore, the values of T n d given in Table | should be considered only as rough estimates. 

The presence of a resonance before a quasi-bound state appears is not surprising. With 
increasing attraction the poles of the S'-matrix should move in the complex plane from the 
resonance area to the quasi-bound state area. This is exactly what our calculations indicate. 
In Ref. || we showed that such transition of the pole happens when Rea^jv changes from 
0.25 fm to lfm. Here we found more exactly that it happens at Re 0,^=0.733 fm. 

The resonant behaviour of rjd elastic scattering should be seen in various processes in- 
volving rjdsjstem in their final states, such as 'yd — > drj and np — > dr\. Recent measurements 



12] of the //-production in the np collisions reveal a bump of the cross-section for the re- 



action np — > drj at a cm. energy below 5MeV. If we suppose that the energy dependence 
of this cross-section is mainly determined by the final state interaction, then this bump 
can be explained by the existence of an rjd resonance. Moreover, since the bump was ob- 
served below 5 MeV, the resonance positions given in Table | imply the rough lower bound 
Re a V N > 0.75 fm. For a reliable estimate, however, one has to perform an explicit calculation 
of the corresponding cross-section. A recent analysis [^J of this reaction with a non-resonant 
final state interaction is consistent with the data of Ref. |L2|| only if Rea^jv = 0.30 fm. Fi- 
nally, it is interesting to note that there is tentative evidence for a similar bump in the 
low-energy region in the reaction ^d — > Xn 
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FIG. 1. Real part of the ?7-deuteron phase-shift as a function of the collision energy. The five 
curves correspond (starting from the lowest one) to Rea V N= 0.55 fm, 0.65 fm, 0.725 fm, 0.75 fm, 
and 0.85 fm. 
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FIG. 2. Imaginary part of the ?7-deuteron phase-shift as a function of the collision energy. The 
curves correspond to the same choices of the Rea V N as in Fig. 1. Within 1 MeV to 6 MeV, where 
the curves do not intersect, the lowest one corresponds to Re a V N=0.85 fm, and every next curve 
above it corresponds to a smaller Rea^jv- 
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FIG. 4. Total cross-section (integrated over the angles) for elastic r/-deuteron scattering as a 
function of collision energy. The five curves correspond (starting from the lowest one) to Rea V N= 
0.55 fm, 0.65 fm, 0.725 fm, 0.75 fm, and 0.85 fm. The dashed line indicates the deuteron break-up 
threshold. 
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FIG. 5. Argand plot for the r]d elastic scattering amplitude in the energy interval from 
to 22MeV. The five curves correspond (from right to left) to Rea V N= 0.55 fm, 0.65 fm, 0.725 fm, 
0.75 fm, and 0.85 fm. When the energy increases the corresponding points move anticlockwise. 
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TABLES 



tiea V N (im) 


E vd ( MeV ) 


T vd (MeV) 


0.55 


8.24 


9.15 


0.65 


7.46 


8.45 


0.675 


7.14 


7.61 


0.70 


6.79 


6.90 


0.725 


6.41 


6.31 


0.75 


6.01 


5.87 


0.85 


4.39 


5.79 


0.90 


3.73 


6.81 



TABLE I. Energy and width of the rjd resonance for various choices of Rea V N. 
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